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Abstract 

We investigate the large scale cosmological perturbation in the Universe with multiple perfect 
fluids. Using the long-wavelength approximation with Hamilton-Jacobi method, we derive the for- 
mula for the gauge invariant comoving curvature perturbation. As an application of our approach, 
we examine the large scale perturbation in a brane cosmology. 
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I. INTRODUCTION 



The analysis of the large-scale cosmological perturbation is important issue to obtain 
the information on the initial density fluctuation which was generated during the era of the 
inflationary expansion of the Universe. The scales of cosmological interest such as large- 
angle cosmic microwave background (CMB) have spent most of their time far outside the 
Hubble radius and have re-entered only relatively recently in the history of the Universe. 

If we concentrate on only the evolution of the large-scale fluctuation, the long- wavelength 
approximation (gradient expansion) of Einstein's equation with the Hamilton- Jacobi (HJ) 
method 0, [3], |J] is convenient to obtain the solution of the fluctuation. Following this 
method, the large-scale fluctuation is given by the lowest order of the spatial gradient ex- 
pansion of the Einstein equation and can be recognized as the fluctuation of the spatially 
dependent integration constants contained in the spatially inhomogeneous "background" 
solution. The large-scale fluctuation of the metric and the matter field can be obtained by 
differentiating the solution of the background Einstein equation with respect to the constant 
of integration contained in it. This approach provides an easier method to give the solution 
of the long-wavelength perturbation and is useful especially for the Universe with multiple 
scalar fields or multiple fluids. 

In this paper, we consider the behavior of the long-wavelength curvature perturbation. 
We aim to derive the formula for the gauge invariant curvature perturbation for the Einstein 
gravity with multiple fluids using the long-wavelength approximation with the HJ method. 
Contrary to the case of multiple scalar fields [0, @, we can obtain explicit scale factor de- 
pendence of the curvature perturbation. As an application of our formalism, we consider 
the large-scale fluctuation in the brane cosmology^, ^ |7|. The analysis of the perturba- 
tion in the brane cosmology is complicated because the perturbation equations do not form 
a closed system on the brane in general. But for the large-scale scalar type fluctuation, 
the perturbation equations contain a closed form on the brane[[5], 01], which may be solved 
without solving for the bulk perturbations. Bulk effects produce a non-adiabatic mode of 
fluctuation caused by the degrees of freedom of the Weyl curvature. The Universe can be 
treated as the system with multi component perfect fluids. 

The plan of this paper is as follows. In Sec. II, we introduce the long-wavelength approx- 
imation of Einstein's equation with perfect fluid. In Sec. Ill, we derive the expression for 
the gauge invariant curvature perturbation using the HJ method. We apply this method 
to the brane cosmology in Sec. IV. Sec. V is devoted to summary. We use units in which 
c = h = 8nG = 1 throughout the paper. 

II. LONG- WAVELENGTH APPROXIMATION OF THE EINSTEIN GRAVITY 
WITH PERFECT FLUIDS 

We consider the Universe with two component perfect fluids of which equation of state 
is given by pW = (i -1 , — l)p^ (i = 1, 2) where F{ is assumed to be constant. The form of the 
line element is 

ds 2 = -(N 2 - WN^dt 2 + 2N i dtdx i + e 2a dx 2 , (1) 

where N is the lapse function, N l is the shift vector and e a is the scale factor of the Universe. 
We assume that each component of the metric can have spatial dependence. By introducing 
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the velocity potential Xii t ne Hamiltonian which derives the long- wavelength dynamics of 
the Universe isffl 



H T = J d 3 x(NH + WHi), 



„— 3a 

H=- e —P* + e 3 "(pV+p(% p« = (e- 3 «P Xl ) A , = (e^P X2 ) r2 (2) 



The three velocity of the fluid is given by 
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e^' is the specific enthalpy of the fluid which has the following relation to the number density 
n (i). 

Using these relations, the specific enthalpy can be written as 

& = r i { P ®f i r, (5) 

and the three velocity is 

By variating this Hamiltonian with respect to a, P a ,Xi, P Xi i N, N l , we have the following 
equations 

5P a 6 3 ' ' v ; 

X, = §^ = Nne*** (P x f^ - Ni Xi>j , (8) 

-Pa = ^ = j e" 3a P Q 2 + 3e 3 « ((A - l)p« + (r a - l)p^) + {N*P a ) d , (9) 

-4 = H = (^Xi)j, (10) 

H = 0, Hi = 0. (11) 

Equations ([11]) are the Hamiltonian constraint and the momentum constraint. If we assume 
that the shift vector is zero, we have 

3(£f = P (1) +P (2 \ (12) 

n / ■ \ ■ / • \ 2 



-JJ^)- 3 {JJ) =(A-i)„« + (r 2 -i),<*>, (13) 

p (1) + 3dAp (1) = 0, p (2) + 3dP 2 p (2) = (14) 
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Although these equations have the same form as for a flat Friedman-Robertson- Walker 
(FRW) universe with perfect fluids, variables a, p^\p^ can have the spatial dependence and 
this system represents the evolution of the Universe with the long- wavelength inhomogeneity. 
The spatial dependence of each variable is determined by the momentum constraint. 

We treat this system using HJ method. By introducing the generating functional 
S[ct, Xii X2], the conjugate momenta of the dynamical variables are replaced by 



P _ SS _5S_ 

a ~5a> Xl ~ 5 X i 



P = ™_ 

X2 5x2 



The Hamiltonian constraint and the momentum constraint become 



H = e 



-3a 



1 f5S 



12 \5a 



+ e 



6a 



-3a 



5S_ 



A 



+ e' 



60 



-3a 



5S 



r 2 



1 /5S\ 



Hi = - — - — a 



3 \5a J 



5S\ 



5a J 



5S 
8X1 



5S 

5X2 



5X2 

X2,i = 0. 



= 0, 



(15) 

(16) 
(17) 



For N l = 0, the evolution equations are 

N 



-6a 



a = e 

6 



Xi — NTi e 



-3a 



5S\ 
5a ) 
5S 



5xi 



A-l 



X2 



-3a 



5S_ 

5X2 



r 2 -i 



By assuming the following form of the generating functional 

S = - 2 J d 3 xe 3a H( X uX2), 



the Hamiltonian constraint becomes 



3H 2 + (-2H Xi yi + (-2H X2 ) r > = 0. 



(18) 
(19) 

(20) 
(21) 



As this equation does not contain xi and X2 explicitly, the form of its solution can be written 

as 



H = H(xi + di,X2 + d 2 ), 



(22) 



where di(x) and d 2 (x) are spatially dependent constants of integration. The momentum 
constraint becomes 

Hi = 2e 3a [-H ti + H xlX i,i + H X2X 2,i] = -2e 3a (H dl d ht + H d2 d 2 ,) = 0. 

For A^* = 0, the conjugate momentum P Xl and P X2 become constant in time and we have 



P 



XI 



-2e 3a H- 



-2 Cl , P 



X.2 



-2e 3a H. 



-2c 2 , 



(23) 



where c±(x) and c 2 (x ) are spatially dependent constants. Hence the momentum constraint 
gives the following relation between spatially dependent constants : 



c\d u + c 2 d 2 i = 0. 



(24) 
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The evolution equations are 

a = NH, 

Xi = iVA(-2if xl ) ri - 1 , (25) 
X 2 = Nr 2 (-2H X2 ) r *-\ 

and ( pip are basic equations to analyze the long-wavelength dynamics of 
the Universe. The spatial inhomogeneity of the dynamical variables a,xi,X2 comes from 
the spatially dependent constants ci(x), C2(x), di(x), d,2{x). Solving HJ equation is formally 
equivalent to obtaining the homogeneous background solution. Once the background solu- 
tion with constants of integration is obtained, the spatial dependence of the long-wavelength 
fluctuation is determined through the momentum constraint. 



(24 



III. LINEAR PERTURBATION ABOUT A FLAT FRW UNIVERSE 

To obtain the curvature perturbation, we consider the linear perturbation about the 
homogeneous background. This means that variables in the equations , (p3|) , ([24]) and 
(p5|) are linearized as follows : 

a -»• a(t) + 8a(t, x), N -> 1 + n(t, x), xi — ► Xi(t) + Sxi{t, x), 
Cj — > Cj + 5ci(x), di — > di + Sd{(x). 

The background equations are given by 

-3H 2 + (-2H Xl ) r ^ + (-2H X2 ) r - = 0, 

a = H, xi = r 1 {-2H xl )*-\ X 2 = r 2 {-2H X2 ) D >-\ (26) 

These equations describes the evolution of spatially flat FRW universe. The constants 
contained in these equations do not have spatial dependence. By using equation (p3|), the 
background solution of the velocity potential can be written as 

Xi = Xi(«> c i> c 2, di, d 2 ), Xi = X2(«, c u c 2 , d ly d 2 ). (27) 
The evolution equations for linear perturbation are 

5a = nH + H Xl Sxi + H X2 5x2, (28) 
*Xi=nA(-2ff Xl ) n - 1 

- 2A(A - l)(-2H Xl ) r ^- 2 (H XlXl 6xi + H X1X2 6 X 2 + H xidl 6di + H Xld2 6d 2 ), (29) 
5 X 2 = nr 2 (-2H X2 ) r >- 1 

- 2r 2 {r 2 - \){-2H X2 ) r *- 2 {H XlX2 5xi + H X2X2 5 X 2 + H X2dl 5di + H X2d2 6d 2 ). (30) 

In equation (f28|), we do not have terms which contain 6di because this equation give the 
relation between 5a and 5xi, and is not the evolution equation. Indeed this equation has 
the same form as the momentum constraint in the conventional cosmological perturbation 
theoryg | . 

We solve the perturbation equation in 5a = gauge. This gauge condition corresponds 
to the zero curvature gauge [0. In this gauge, by eliminating the lapse function n and using 
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the background solution 
form: 



we can rewrite the evolution equation for 5\i as the following 



5xi + Sdi 
6x2 + 6d 2 



X n X 



i (Sxi + 6dA , H xi 8di + H X2 5d 2 /r x (-2F xl ) r i 



5x2 + 6d 2 



+ 



H 2 



r,(-2H X2 ) r ^ 



(31) 



where the matrix X is defined by using the background solution ( p7p as follows 



Xl,Cl Xl,C2 

X2, C1 X2,c 2 



(32) 



If we substitute the following form of the solution to the evolution equation (|31"D , 



+ 8di 
5x2 + 6d 2 



X 



we have 



X 



Ci 

c 2 



-3a 



c - 



H 



Ci(t) 

c 2 (t) 



Xl,a 
X2,a 



(33) 



(34) 



where c = C\5di + c 2 5d 2 is a constant which has no spatial dependence. This comes from 
equation fl24|). Therefore, an d <^X2 are given by 



6xi 
6X2 



X 



5c\ 
5c 2 



+ Co da 



H 



X 



1 / Xl,c 
X2,c 



8d x 
5d 2 



(35) 



where 5c\(x) and 5c 2 (x) are spatially dependent constants. 

Now we evaluate the gauge invariant variable corresponding to the curvature perturbation 
on comoving slice 



K 



-5a — H 



(p(i) +p (i)yi) + (p(2) +p (2)y2) 

(p(l)+p(l)) + (p(2) + p(2)) 



-5a + -r-(H Xl 5xi + H X2 5x 2 ), 



(36) 



where we have used equation (6) to express the velocity perturbation by 5\- This quantity 
gives the spatial curvature perturbation on the comoving slice and equivalent to Bardeen's 
parameter ( in the long- wavelength limit. To proceed our calculation, we prepare formulas 
which comes from the conservation law equation (23). By differentiating equation (^) with 
respect to a, we have 



M 



Xl,c 
X2,c 



H X1 
H X2 



M 



^XlX2 H-X2X2 



(37) 



By differentiating equation (|23|) with respect to c\,c 2 , we have 

MX o-< ( J J 
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Using these relations, the gauge invariant curvature perturbation 1Z in the 5a = gauge 
becomes 



n=S-(H xl 5 X i + H X2 6 X 2) 

ti 

where energy density of each fluids is given by p^ = (— 2c 1 ) ri e~ 3riCt , p^ = (— 2c 2 ) r2 e~ 3F2a . 
This is the main result of this paper. The first term of this expression corresponds to the 
contribution of the growing mode of adiabatic and iso-curvature perturbation. This term 



is same as the one derived by using the local conservation of energy momentum tensor |p 
The second term corresponds to decaying mode. In our analysis, the decaying mode cannot 
have spatial dependence, but it gives the correct time dependence. For A = A , the system 
reduces to a single component fluid case and 7Z obeys the following evolution equation 

K = for A = 1, (39) 

It + 3aTZ = for A ^ 1 (40) 

In the situation when the energy density of one fluid dominates p^ 3> p^ 2 \ the time de- 

2 

pendence of the scale factor is e a ~ (t/t ) 3r i and the behavior of the curvature perturbation 
is given by 

6a 1 (5c 2 5a\ r 2P2 ( f ^_ 3a H^_ 3a 



a 3 V c 2 a J Api \J H 
5a 1 ( 5c 2 5a\ A(-2c 2 ) ra w „_ 3(A _ A)a 2(7 1 - Vj _ 



, x e -^'^ a - c t x v 1 ^ e W^)«. (41) 

a s\c 2 a J A(-2cO r i 2- a 

The first term corresponds to the adiabatic growing mode which is constant in time. The 
second term corresponds to the iso-curvature mode and the third term is the decaying mode. 
The decaying mode vanishes in the A = 1 case. 



IV. APPLICATION TO THE BRANE COSMOLOGY 

We apply the method developed in Sec. Ill to evaluate the large scale curvature perturba- 
tion in the brane cosmology 0, [| [7|. By assuming the Randall- Sundrum type brane model, 
the four dimensional Einstein equation on the spatially homogeneous brane universe is 

SH 2 = p+^p 2 + p £ , 

-2H-3H 2 =p+ ^(p 2 + 2pp) + 2L % (42) 
p + 3H(p + p) = 0, p = (r-l)p, p £ + 4Hp £ = 0, 

where I 2 is the inverse of brane tension and p £ is the energy density of the Weyl matter. 
The unconventional feature of the brane cosmology is the appearance of p 2 and pe terms 
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on the right hand side of the Einstein equation. The Hamiltonian for the long-wavelength 
dynamics which derives equation in the homogeneous limit is 



H T = J d 6 x(NH + N l Hi), (43) 
^ = -^ 2 + e 3Q (p+^p 2 + P,), p=(e- 3 «P Xl ) r , p £ =(e^P X2 )K (44) 

Pai 

Hi = -y- - p a&,i ~ P X iXx,i - P x2X2,t- (45) 
The effective energy density and the effective pressure for the perfect fluid are defined by 

I 2 

Pm= P + y^P 2 , (46) 

PM = (r-l)p+^(2r-l)p 2 . (47) 
The specific enthalpy for the effective matter field is given by 

r-i / I 2 2 

t M = Fp— ( 1+ gP 

For the Weyl matter, 

4 1/4 

Hence the three velocity of the effective matter field and the Weyl matter becomes 

8xi s Xi 



vm = ~ 



V£ 



es rp^ (1 + f p) 

^X2 8x2 



The gauge invariant variable corresponding to the curvature perturbation on comoving slice 
is 

(Pm + Pm)vm + (Pe + Ps)v£ 



n=-5a-H 
= -5a + H 



(Pm+Pm) + (p£+Ps) 
p^dxi + p£ 3/ *$X2 
r(p + f P 2 ) + iPe 



= - 6a + ^( Hl S Xl + H 2 5x2). (48) 
H 

This expression is same as the result of standard cosmology fl3"6|). The explicit form obtained 
in the 5a = gauge is 

* - A^4M^ + J[ *- - M ■ <-o 



3 r(p + ^) + l P£ U V7 h 

The first term of this expression reproduces the result of Ref. || which was derived by using 
the local conservation law of the fluids. Our method include the contribution of the decaying 
mode. The time dependence of the curvature perturbation can be obtained in the following 
three cases: 
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a. l 2 p 2 dominates The time dependence of the scale factor is e a ~ (t/to) 3 ^ an d 
6. p£ dominates The time dependence of the scale factor is e a « (t/t ) 1//2 and 



<5c 2 






5ci 


3c 2 




V c 2 


Ci 


5c 2 




^c 2 




3c 2 




V c 2 


Cl 



(~ 2Cl ) r x p (4-3r)a , 12 (- 2 Cl) 2r -2(3r-2)a\ /ci \ 

,(-2p a )Vs Xe + 6(-2c 2 )V3 Xe J ^ 

- c t x e~ Q . 



2 

z. p dominates The time dependence of the scale factor is e a ~ (t/t )zr and 



<5ci 


- 4 i 


^c 2 


8ci 


3ci 


9r 


V c 2 


Cl 


5ci 


. 4 1 


^5c 2 




3ci 


9r 1 


V C 2 


Cl 



/a 



|^ x e-^n _ Coto x ^_l) e -I(2-r)« (52) 



V. SUMMARY 



We have derived a formula of long-wavelength curvature perturbation for the Universe 
with multi-component fluids using HJ method. The non-decaying mode of the curvature 
perturbation is explicitly expressed as the function of scale factor of the Universe. This mode 
is the mixture of the adiabatic growing mode which is constant in time and the time varying 
iso-curvature mode. In the case of the multiple-scalar field]!]], the curvature perturbation 
cannot be expressed explicitly as the function of the scale factor (see Appendix) because we 
cannot obtain the solution of the background scalar fields as the function of the scale factor 
in general. For the perfect fluids case, each component of fluids evolves separately and this 
makes it possible to obtain the background and the perturbation solution. 

Although the decaying mode obtained in this paper reproduces the same time depen- 
dence as the result of the conventional cosmological perturbation theory, it cannot have the 
spatial dependence because of the momentum constraint. The homogeneous decaying mode 
obtained here can be absorbed to the background solution by re-defining the background 
constant of integration. This is related to the feature of decaying mode in zero curva- 
ture gauge and co-moving gauge ; the decaying mode under these gauge conditions becomes 
higher order in long- wavelength expansion compared with that of other gauge conditions [|l2j . 

As the application of the HJ method, the large scale curvature perturbation in the brane 
cosmology was derived. Our method reproduced the result obtained by using standard cos- 
mological perturbation theory with the local conservation of the energy- momentum tensor 0. 
It is interesting to investigate the brane model with the scalar field. But we could not obtain 
the Hamiltonian which derives the long-wavelength dynamics of the brane cosmology with 
the scalar field. This is left as our future problem. 
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APPENDIX A: CURVATURE PERTURBATION IN THE MULTIPLE SCALAR 
FIELDS SYSTEM 



We review the derivation of the long-wavelength curvature perturbation in the Einstein 
gravity with multiple scalar fields 0]. HJ equation and the evolution equation are 

- 3H 2 + 2Hl + 2Hl + Vfa, <h) = 0, (Al) 
a = NH, (A2) 
<j> 1 = -2NH th , fa = -2NH <h (A3) 

The solution of the HJ equation is given by H = H(<px, <p 2 , d±, d 2 ) where di(x) and d 2 (x) 
are spatially dependent constants. From the momentum constraint, the spatial dependence 
of these constants are determined: 

Hfadij + H d2 d 2 A = 0. 

Other constants of motion in the system are given by 

e 3a H dl = ci(x) = e - 3ao ^, e 3a H d2 = c 2 (x) = e" 3 ^/^), (A4) 

From this, we can express 0i, <p 2 as the function of the scale factor: 

0i = (j)x(a + a ,f, d 1} d 2 ), (j) 2 = (j) 2 (a + a , f, d h d 2 ). (A5) 

The evolution equation for the perturbation is 

5a = nH + H^Sfa + H^Sfc, (A6) 
54>i = -2nH^ x - 2(H Ml 54> 1 + H^Jfa + H^Sdx + H^ ld2 5d 2 ) : (A7) 
502 = —2nH^ - 2(H^ 2 5(j)i + H^^Sfe + H^Sdi + H t f >2d2 5d 2 ). (A8) 

The gauge invariant variable corresponding to the curvature perturbation on comoving slice 
is given by 

H 3 

1Z = —5a — 7777- (0i, a 50! + (t> 2 ,J<p 2 ). ( A9) 

lti a 

We evaluate this expression in the 5a = gauge. For this purpose, we prepare the relation 
derived form equation By differentiating equation (|A4|) with respect to c\, c 2 , we have 



MX = e~ 3a l, M = ( ) , X = ( l]' Cl ™ ) (AlO) 



H^ ldl H t 

The evolution equation in the 5a = gauge can be written as 



V n " WJ ~ H W dl H 2d2 J Ud 2 ) • 



By substituting 



Cx(t) 
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to the evolution equation, we have 



X 



Co 



H 



M 



T 



5di 
5d 2 



2c 



-3a 



H 



(X 



-1\T 



5d\ 
5d 2 



and the solution of equation ( |A1 1| ) is given by 



Hi 

5<j)2 



X 



Sci 
5c 2 



-3a 



+ / da- 



H 



[x-^x- 



-1\T 



Sdi 
5d 2 



(A12) 



where 5c\,5c 2 are spatially dependent constants. Using this solution, the curvature pertur- 
bation becomes 



K 



H 3 



2H n 



5c\ 
Sc 2 



-3a 



For the growing mode 5d 



l,a (j>2,a) X { , ! j : 

5d 2 = 0, we have 
7Z = Sao 



da- 



H 



-1\T 



5di 
5d 2 



l,a<Pl,f 



+ 



>2,a92J 



+ 



5f- 



(A13) 



(A14) 



J 2,a, 



The first term corresponds to the adiabatic growing mode which is constant in time, and the 
second term corresponds to the iso-curvature mode. The explicit scale factor dependence 
can be obtained with the assumption of the slow rolling approximation and the specific form 
of the scalar field potential KJ. 
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